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Frequency light shifts resulting from the localization effects and effects of the quantization of trans- 
lational atomic motion in an optical lattice is studied for a forbidden optical transition J=0— * J—0. 
In the Lamb-Dicke regime this shift is proportional to the square root from the lattice field intensity. 
With allowance made for magneto-dipole and quadrupole transitions, the shift does not vanish at 
the magic wavelength, at which the linear in intensity shift is absent. Preliminary estimates show 
that this shift can have a principal significance for the lattice-based atomic clocks with accuracy of 
order of 10" 16 -10" 18 . Apart from this, we find that the numerical value of the magic frequency de- 
pends on the concrete configuration of the lattice field and it can vary within the limits 1-100 MHz 
(depending on element) as one passes from one field configuration to another. Thus, theoretical 
and experimental investigations of contributions originated from magneto-dipole and quadrupole 
transitions are of principal self-dependent interest. 

PACS numbers: 42.50.Gy, 42.62.Fi, 42.62.Eh 



The last few years were marked by a principal ideolog- 
ical [l[ and experimental @, H, 0, la, 13] breakthrough in 
the field of fundamental laser metrology based on a possi- 
bility to work with strongly forbidden optical transitions, 
using large number of neutral atoms confined to an opti- 
cal lattice at the magic wavelength. This great advance 
expressed itself in a transition to spectral structures with 
Hertz 0, @| and potentially sub-Hertz widths, that ap- 
proaches a prospect of the creation of atomic optical 
frequency standards with unprecedented fractional fre- 
quency uncertainty and accuracy at a level 10 -17 -10 -18 . 

However the achievement of such extraordinary metro- 
logical characteristics in real frequency standards is a 
challenging goal. On the way to this goal it will be 
necessary to determine sources and values of systematic 
errors. Therefore now there comes a period, when fre- 
quency shifts of various origin should be thoroughly in- 
vestigated. For example, the second-order (in intensity) 
light shift due to the atomic hypcrpolarizability has been 
studied experimentally [B|, 0] and theoretically [U . 

In the present paper, taking into account contributions 
due to magneto-dipole and quadrupole transitions, we 
investigate a previously unknown frequency light shift 
caused by the quantization of atomic translational mo- 
tion in an optical lattice for a forbidden optical transition 
J g =0— *J e =0 (for instance, 1 So— > 3 Pq in alkaline-earth- 
like atoms) . This shift is proportional to the square root 
(in the Lamb-Dicke regime) from the lattice field inten- 
sity, and it does not vanish at the magic wavelength A m , 



at which the first-order (in intensity) light shift cancels. 
Preliminary estimates show that this unavoidable shift 
has a principal significance for the definition of metro- 
logical characteristics of lattice-based atomic clocks. 

Note, the influence of the electric quadrupole (E2) and 
magnetic dipolc (A/1) effects on the clock levels in Sr 
atoms was first discussed in [l|. The calculations have 
demonstrated that the E2 and M 1 dynamic polarizabili- 
ties at the magic wavelength is of order 10 -7 of the elec- 
tric dipole (El) polarizability. However, the calculations 
in [l[ referred to the usual plane- wave approach which did 
not account for the specific inhomogeneous distribution 
in space of electric and magnetic field components of the 
optical lattice. The effects of inhomogeneous field distri- 
bution of a standing wave are considered in the present 
paper in detail. 

Consider an atom confined to an optical lattice, which 
is induced by a one-dimensional elliptically polarized 
standing wave (with the frequency uj). The electric field 
vector has the form: 



E(r, t) = E e cos(k • r) 



c.c. 



(1) 



where Eq is the scalar amplitude; k is the wavevector 
(k=\k\—uj/c); e is the complex unit polarization vector 
(e • e*)=l. The condition (e • k)=0 is satisfied due to the 
transversality of electromagnetic field. 

First we consider the frequency shift of a transition 
J g =0— >J e ~Q in a potential induced only by the contri- 
butions of electro-dipole transitions Jj=0— >J=1 (j=g,e). 
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For the standing-wave field ([T]) the light shift (potential) 
of a j-th level is spatially modulated and it has the fol- 
lowing form: 

Uf{v) =-W jC os 2 (fcz); W,.>0 (j = g,e). (2) 

Here, for convenience, we direct the z axis along the 
wavevector k. The potential amplitude W» depends on 
the frequency u and it is proportional to |i?o| 2 - 

Then we quantize the translational motion of atoms. 
To do this we will assume the Lamb-Dicke regime, 
when atoms are localized in the field antinodes kz=lir 
(7=0, ±1, ±2...) on the size much less than the wave- 
length. In this case we can describe the atomic motion 
in the harmonic approximation. For the sake of dcfinitc- 
ness, we will work near the point 2=0. At the condition 
|fcz|<Cl we can use the approximation cos 2 (/cz)«l— k 2 z 2 , 
which allows us to write the shift ^ as the harmonic 
oscillator potential: 

„f to „- W , + (,_,,.), (3) 

where M is the atomic mass, and the oscillator frequency 
f2j for j-th level has the form: 

1 / 2W ■ k 2 
^• = ^\/^-; 0=*e)- (4) 

The following relationships are obvious: 

WjOtl, oc v7, (5) 

where I—c\Eq\ 2 /2tt is the field intensity in the lattice 
antinode. 

Using for the potential ([3]) the standard theory of har- 
monic oscillator, we write the energies of the upper and 
lower levels with consideration for the vibrational struc- 
ture: 

£ s (n) = £f ] - Wj + Ml, (n + 1/2) , (j =g,e), (6) 

where h—2irh; is the energy of unperturbed j-th 

state in a free space; n=0,l,2,... is the vibrational quan- 
tum number (see in Fig.l). 

Consider the optical frequency of transition between 
vibrational levels with the same quantum numbers n: 



i/W - (W e - W g )/h + (fi e - Q g )(n + 1/2) , 

where u^ — {£^—£^)/h is the frequency of unper- 
turbed Jg—0— ► J e =0 transition. Taking into account the 
relationships l|5|), the frequency shift can be written as: 

Av nn = v nn - = a(u)I +(n+ l/2)P(w)Vl , (8) 



n 




FIG. 1: An illustration of optical transitions between vibra- 
tional levels for atoms confined to an optical lattice. 

where the coefficients a(u>) and /3(u>) are individual char- 
acteristics of chosen element. They are defined as follows: 

a(uj)I = -(W e -W g )/h; (3(uj)Vl = n e - n g . (9) 

Thus, despite of the fact that we consider the potential 
([2|), which is proportional to I, the effects of the quanti- 
zation of atomic motion lead (in the Lamb-Dicke regime) 
to the appearance of the additional square root depen- 
dence ((xVl) for the frequency shift ([8]). Note that be- 
yond the Lamb-Dicke regime the intensity dependence of 
the frequency shift is more complicated due to the an- 
harmonicity. 

In the case of potential ([2|), which is induced by 
electro-dipole transitions only, at the magic frequency 
w m =2-Kc/\ m we have simultaneously a(u> m )—0 and 
(3{bj m )=Q, because W e (u} m )=W g (uj m ). However, if we 
take into account contributions due to magneto-dipole 
and quadrupole transitions, it will not be so, i.e. 
/3(wm)7^0, that can have a principal importance for 
metrological characteristics. Let us prove this. 

For the standing-wave field {1} the magnetic field vec- 
tor B has the form: 

B(r, t) = B Q e B sin(k • r) e^* + c.c. , (10) 

where B =iE is the scalar amplitude of the magnetic 
field, and e B = [kxe]/fc is the unit polarization vector 
of the magnetic field. As a result, the contribution to 
the potential of j-th level (oc|B| 2 ) due to magneto-dipole 
transitions Jj=0— >J=1 has the spatial dependence dif- 
ferent from the potential |2]): 

Uf{v)=B j sin 2 (kz); (j = 9,e), (11) 
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where the potential amplitude Bj is proportional to the 
intensity /, but the frequency dependence (on u>) of Bj 
differs from It can be easily shown that the con- 

tribution due to quadrupole transitions Jj=0^J=2 has 
the same spatial dependence (for ID standing wave): 



Uf(r) = sin 2 (fa 



(j =9,e). 



(12) 



Hence, for j-th level the total potential proportional to 
the intensity / has the form: 



U j {T) = Uf{v) + Uf{v)+uf{v) 



(13) 



-Wj cos 2 (fa) + {B } + Qj} sin 2 (fa); (j = g, e). 

As is seen, the spatial dependencies of potentials U g {r) 
and Z7 e (r) for the lower and upper levels of the forbidden 
J g =0— > J e =0 transition are different due to the contribu- 
tions induced by magneto-dipole and quadrupole transi- 
tions, i.e. U g (r)/U e (r)^ const. As a consequence, strictly 
speaking, the condition U e (r)=U g (r) does not hold for 
any frequency. Hence, the ideal notion of the magic fre- 
quency u> m does exist only for a single running wave, 
for which, however, the confining optical lattice poten- 
tial is absent (spatially uniform light shift) and there- 
fore this case is not important for lattice-based atomic 
clock. However, since the contribution due to electro- 
dipole transitions dominates, the other contributions can 
be considered as very small perturbation, which, never- 
theless, can influence the metrological characteristics of 
frequency standards. 

Expanding the expression (| 13[) in (fa) powers and us- 
ing the harmonic approximation (i.e. cos 2 (fa)«l— k 2 z 2 
and sin 2 (fa)«fc 2 z 2 ), we obtain the expression for the po- 
tential analogous to J3]): 



UAr) 



-w, 



M(2tt^) 2 z 2 



(j = 9,e), (14) 



from which all the following formulae ©-© are deduced. 
However instead of ttj (see eq.Q) now we have to use 
the other expression for the vibrational frequency flj: 



J 2tt 



1 /2{W,+S, + Q,}fe 2 



M 



(j = 9,e). (15) 



From the equations (TH| and (fl3|) it follows that for ID 
standing wave the magneto-dipole and quadrupole tran- 
sitions affect only the coefficient 0{u>) in the formula for 
the shifts {8]), while the coefficient a(oj) is governed by 
the electro-dipole transitions solely as before. 

The magic frequency of lattice field Lo m can now 
be determined from the condition of cancelation of 
the linear (with respect to the intensity /) shift in 
([5]) (i.e. a(u> m )=0). In this case we find that 
W e (uj m )=W g (uj m )=W. Then the remaining shift (cx\/7) 
in eq.© differs from zero: 



(n + l/2)/% ro )V7^0. 



(16) 



Expanding the expression (|T5|) in small parameter 
\(Bj+Qj) /W|<Cl and leaving only the first-order term, 
we obtain: 



/?(u; m )Vl = (Q e - O s ) « ft ( V 
Here the frequency f^ -* is equal to: 



2tt V M 



(17) 



(18) 



and its value coincides with the vibrational frequency, 
which is practically the same (at ui m ) for the upper and 
lower levels of the clock transition J g =0^J e —Q. The 
dimensionless small coefficient ry in (I17|) does not depend 
on the intensity / and polarization e, and it is defined as: 



Vb 2W 



(19) 



Q 



9 . 



2W 



\Vb>Vq\ < !■ 



The first term rj B in (| 19|1 is governed by the magneto- 
dipole transitions, and the second term tjq is caused by 
the quadrupole transitions. 

Now let us estimate the metrological significance of 
the unavoidable square-root shift (fT?"]) . Basing on the 
very general estimations, we can expect that for different 
elements the value of the coefficient rj is of order of 1CP 7 - 
1CP 6 . For typical optimal (with respect to the intensity 
/) experimental conditions the vibrational frequency 
is of order of several tens kHz. Then even for ?7^10~ 7 
the shift (fT6|) - (fl~7| is estimated as Af„ n ~0.01 Hz, what 
is principal for frequency standards with the accuracy 
10 -17 -10~ 18 . It is very important that this shift can not 
be substantially reduced by the decrease of the field in- 
tensity in view of weak square-root dependence \fl (in 
contrast with the hyperpolarizability leading to the shift 
oc/ 2 ). Note also that the discussed shift should be taken 
into consideration in experiments on the precision mea- 
surement of the magic frequency uj m . 

The results obtained above can be generalized to the 
case of arbitrary field configuration (including 2D and 
3D optical lattices), when the electric field vector has 
the general form: 



E(r, t) = E(r)e-' iwt + c.c; E(r) = ^ E a e 



i(k„r) 



(20) 



where E a is the vector amplitude of a-th running wave 
with the wavevector k a (|k a |=fc=a;/c). The spatial de- 
pendence of potential induced by the electrodipole tran- 
sitions is governed by the expression: 

Uf(r)= Wj (u,)\E(r)\ 2 ; (j = g,e), (21) 

where the frequency dependence wAuj) is an individual 
characteristics of chosen element. 
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The potential induced by the magneto-dipole transi- 
tions has the form: 

Uf(r)=b j (cj)\B(r)\\ (j = g,e); (22) 

B(r) = J2 e l(k " r) [n ka x EJ , (n ka = k a /k) . 

a 

And, eventually, the contribution due to quadrupolc 
transitions can be presented in the form of the follow- 
ing scalar product: 

(r)=r»(Q a .Q$); (j = g,e), (23) 

where the covariant components of irreducible tensor of 
the second rank Q 2 are written as: 

Q 2 , = I> i ( k ° r >{n ka ®E a } 2 ,, (g = 0,±l,±2). (24) 

a 

The definition of the tensor product of two arbitrary vec- 
tors {&(S>h} 2q , and the expression for the scalar product 
of the type (|a£§)b} 2 -{c£§)d} 2 ) can be found, for example, 
in the book [9|]. 

Evidently, that in the general case all the spatial de- 
pendencies (|2"Tj) . (|22p and (|2"3"|) are different from each 
other. In this case the magneto-dipole and quadrupole 
transitions can also affect the linear term in ([8]), i.e. the 
coefficient a(oj). This will take place, if in minima points 
{r m m} of the electro-dipole potential U^(r) we have 
{fmin)^ and/or (r min )^0 (as distinct from the 
ideal one-dimensional standing wave ([1])). As the result 
the magic frequency uj m will differ from its value for the 
ID standing wave, for the later it is worth to introduce 
a separate notation uj^. In the order of magnitude the 
region of possible variations of uj m for different field con- 
figurations can be roughly estimated as ^r/uj^, i.e. for 
some elements one can expect the variation at a level 100 
MHz (or even more). Even in the one-dimensional case, 
when the counterpropagating waves are unbalanced (i.e. 
when they have different amplitudes) the magneto-dipole 
and quadrupole contributions have influence on the coef- 
ficient a(u>) in ([8]), changing the magic frequency u> m . 

Concluding, in the present paper for a strongly forbid- 
den optical transition J=0— >J=0, taking into account 
contributions due to the magneto-dipole and quadrupole 
transitions, we have investigated the previously unknown 
frequency light shift caused by the quantization of atomic 
translational motion in an optical lattice. Using ID 
standing wave as an example, we have proved that the 
above mentioned factors leads to the existence of the ad- 
ditional frequency shift, which has square-root depen- 
dence on the lattice field intensity in the Lamb-Dicke 
regime. It has been also shown that this shift does not 
vanish at the magic wavelength. The main intrigue is 
that the new square-root frequency shift can have a prin- 
cipal significance for the metrological characteristics of 



atomic clocks and, due to this reason, it should be thor- 
oughly investigated. In view of the weak square-root de- 
pendence ocvT and small absolute value, it maybe dif- 
ficult to measure experimentally this shift with a good 
accuracy (especially against the background of strong de- 
pendence oc I 2 due to the hyperpolarizability), what in- 
creases the importance of theoretical calculations. Apart 
from this, it has been found that the numerical value of 
the magic frequency depends on the concrete configura- 
tion of the lattice field and it can vary within the limits 
I-I00 MHz (depending on element) as one passes from 
one field configuration to another. 

The obtained results can be significant for the choice 
of experimental methods in the case, when the unavoid- 
able shift turns out to be important. For example, from 
eq. tfTB"]) it follows that the shift Ai/ nn depends on the vi- 
brational quantum number n. Therefore, before the be- 
ginning of spectroscopic measurements it is advisable to 
transfer atoms to the lowest vibrational level with n=0 
(similarly to the method successfully realized in [lfj). Be- 
sides, it is preferable to use a resonator for the formation 
of ID standing wave for the purposes of more reliable 
balance of counterpropagating waves. 
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